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Abstract 

Given a set S of points in the plane, a geometric network for S is a graph G with vertex set S 
and straight edges. We consider a broadcasting situation, where one point r £ S is a designated 
£SJ source. Given a dilation factor 5, we ask for a geometric network G such that for every point v £ S 

t-H there is a path from r to v in G of length at most <$|ru|, and such that the total edge length is 

minimized. We show that finding such a network of minimum total edge length is NP-hard, and give 
an approximation algorithm. 
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1 Introduction 



Given a set S of points in the plane, a geometric network for S is an edge-weighted graph G with vertex 
set S and straight edges. The weight of an edge (u, v) is the length of the segment uv, that is, the 
Euclidean distance \uv\ of the two points. 

Various types of networks, such as communication networks, road networks, or telephone networks, 
have been modeled as geometric networks. One important parameter of a geometric network is its total 
cost t{G): the sum of all edge lengths. The network that minimizes the cost while connecting all points 
in S is the Euclidean minimum spanning tree of S. Another well-studied parameter is the dilation of a 
network. For two points u,v £ S, the dilation of the pair (it, v) is defined to be the ratio 

cn 

^ a , \ d G (u,v) 

\uv\ 

V.Q of the length of the shortest path between u and v in G and the distance \uv\. The dilation of a network 

is commonly defined as the maximum of A.q(u,v) over all pairs u,v G S. The network minimizing the 
dilation is the complete graph, which has dilation 1. 

The complete graph has prohibitively large cost, while the minimum spanning tree may have large 
dilation. Balancing these two parameters has been the subject of much research in the literature, we 
refer to the book by Narasimhan and Smid [4] for an overview. 

In this paper, we consider a broadcasting situation, where one point r € S is a designated source, and 
the purpose of the network is to broadcast information from r to all the other nodes. A node v receives 
the information with delay dc(r, w), and we are interested in the relative delay, which is the dilation 
A(j(r, v) = dc{r, v)/\rv\. We define the (relative) delay of G to be the maximum 

\tn\ a ( \ d G (r,v) 

A G = max An(r,v) — max — , — , — . 
»ES\{r} «eS\W \rv\ 

Note that we can assume our network G to be a tree, as the shortest-path tree with source at r will have 
the same delay as G itself. 
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The network minimizing A(G) is the star consisting of all edges (r, v), for v € 5\{r}. It has A(G) = 1, 
but its cost £(G) may be prohibitively large. The minimum spanning tree, which minimizes £{G), may 
have large delay. In practice, a trade-off between both factors needs to be achieved. Related notions have 
been studied in the literature. For instance, computing the tree T that minimizes £{T) while bounding 
the distance of all nodes from the source in the graph is NP-complete [5j. Similarly, minimizing £{T) 
while keeping all nodes at most k hops from the source is NP-complete pQ. 

Our problem, which we term the single-source dilation-bounded minimum spanning tree problem, 
takes as input the point set S, the designated source point r € S, and a delay bound 5^1. We ask for 
a spanning tree T of smallest possible cost £(T) such that A(T) ^ 5, or, in other words, such that for 
every v <E S \ {r}, we have dx(r, v) ^ S\rv\. 

When we set 6=1, then the answer is simply the star connecting r with all other points. When 
we set 5 ^ n — 1, then the answer is the minimum spanning tree of S, and so our problem interpolates 
between these two networks. 

We show that solving the problem exactly is NP-hard by reduction from Knapsack in Section [3j This 
leads us to consider an approximation algorithm. Since this algorithm is rather simple, we present it 
first in Section [2j 

2 Approximation Algorithm 

An approximate solution to our problem can be computed using known algorithms for <5-spanners. We 
introduce these results first. 

Given a point set S in the plane, a <5-spanner for S is a geometric network on S such that dc(u, v) 
8\uv\ for each pair of points u, v € S, where 5^1. Many algorithms to compute a <5-spanner for a 
given point set have been given in the literature, see again the book by Narasimhan and Smid [1] for an 
overview. We will make use of a result by Gudmundsson, Levcopoulos, and Narasimhan [3J, who showed 
the following: 

For any real number 6 > 1, a <5-spanner of a point set S in M. d can be constructed in O(nlogn) time 
such that the spanner has C\n edges, maximum degree C2, and cost c$i{M), where M is the minimum 
spanning tree of S, and ci,C2,C3 are constants depending on S and the dimension d. 

We can now explain our simple approximation algorithm: 

Theorem 1. Given a set of n points S in the plane, a designated source r € S , and a real constant 
5 > 1, we can construct in time 0(n\ogn) a tree T with vertex set S such that the delay A(T) ^ 5, and 
£(T) ^ c$£(M), where M is the minimum spanning tree of S and cs is a constant depending on 5. 

Proof. We first run the algorithm by Gudmundsson et al. [3] to obtain a (5-spanner G for S, with total 
cost £(G) < cg£(M), where eg is a constant that depends on S. 

We then compute the shortest path tree T with source r in G. Clearly we have £(T) ^ £{G) ^ cg£(M). 
For any v € 5\{r}, we have dr(r, v) = dc(r, v) ^ 6\rv\ since G is a (5-spanner, and so the delay A(T) ^ 5. 

G can be computed in time 0(n log n) [3], and the shortest-path tree T can be computed in time 
0(n log n), for instance using Dijkstra's algorithm, using the fact that G has 0(n) edges. □ 

3 Single-Source Dilation-Bounded Minimum Spanning Tree is NP-hard 

In this section, we show that the decision version of our problem is NP-hard. 

Theorem 2. Given a set S of points in the plane, a designated source r € S, a delay bound 8^1, and 
a cost bound K 0, it is NP-hard to decide if there exists a tree T for S with delay A(T) SC 6 and cost 
£{T) < K. 

The proof is by reduction from Knapsack, which we define first. Knapsack is well known to be 
NP-complete [2]. 

Knapsack: Given a set of n items, each having an integer profit pi > and an integer 
weight Wi > 0, as well as a profit bound P > and a weight bound W > 0. Does there exist 
a subset of items with total weight at most W and total profit at least PI 
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Our reduction takes as input an instance of Knapsack, and produces an instance of single-source 
dilation-bounded minimum spanning tree. For a set of n items with profits pi and weight to,, we 
construct a set S of 3n + 4 points. For each item i, we construct three points <ii,&i,Cj. In addition, 
we create three extra points do, d±, c?2, as well as a source r. 
We first define 

a-i ■= Pi + w 4 > 

Pi := 2pi + Wi — Ui + Pi > on 

7, := 3pi +Wi= ^ +pi> fa. 

We thus have 

Q. + Pi = 3pi + 2wi = ji + Wi > ji, 

and so < cti < /3i < ji < ai + /3i. This implies that the three values c^, Pi and 7* satisfy the triangle 
inequality. We also set 

n 

m := max 7^ L := Y^(7i + m), 

l<i<n * — ' 
i=l 

Construction of S. Our construction starts by placing the source r at the origin. We then place the 
2n points a.; and hi on the negative y-axis as follows (see Fig. [T]) : 

a, := (0,-4L-£;~i(7 J+ ™)) 
&i := a>i - (0,7i) 



10L. 



d 2 £ 6L_ 



«i 

Cl 



</() 



3L 



Figure 1: The construction of S 
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We have |a.j&j| = ji and |&jdj+i| = m. We can now place c,, for 1 ^ i ^ n on the right side of the 
y-axis such that \aiCi\ = /?, and \cibi\ = on. Finally, we define do, d\, and d 2 as follows: 



d 
di 
d 2 



(0,-5L) 
(0,-8L) 
(-6L,-8L) 



Regular trees. We classify the edges connecting points of S into regular and irregular edges. Regular 
edges are the edges rai, b n do, dodi, d\d 2 , the edges afii, 6jCj, aiCi, for 1 ^ i ^ n, and the edges ^a^+i for 
1 < i < n. A regular tree is a spanning tree on 5 which contains only regular edges. A tree containing 
irregular edges is irregular. 

In a regular tree T, the dilation of the pair (r, d 2 ) is much larger than the dilation of any other 
pair (r, v), and so we have 

Lemma 3. The delay of a regular tree T is A(T) = A T (r, d 2 ). 

Proof. Let Li — \a\Cbi\ = X]j=i(7i + m )- Since T is a tree, there is a unique path from hi to aj + i and we 
have dr{ai 7 hi) < |ajCj| + \cibi\. Since Qj < /3, < 7* ^ m, we have 

« i— 1 

d T (r, h) < |r, oi| + ^(K C J'I + l c Al) + \ b 3 a 3+i\ 

3 = 1 3 = 1 

i i — 1 

= 4L + ^2(/3j + aj) + ^ m < 4L + L + L i = 5L + L i 

3=1 3=1 
i-1 

dr(r,Ci) < |r,oi| + y^(Kcj| + |c,-6j| + |6j-a i+ i |) + |a<6i| + |bjCi| 

i — 1 2 

= 4L + ^(/3j + to) + aj + 7,) < 4L + Lj + L = 5L + L, 

3=1 3=1 

Note that h) — drir, at) + dr(ai, hi) > dx{r, aj). Thus, we get 

A , \ d T (r,ai) d T (r,bi) 5L + U 

A T (r,a, = — — < — — < < 1.25 

|r-Oj| | raj | 4L + Li 

A , , n d T (r,bj) d T (r,bj) 

Since |rcj| > \rai\, we also have 

A T (r, Ci ) = ^<^ = ^<1.25 
\rci\ \rai\ 4L + L { 

d T {r,d ) d T (r,b n )+m 5L + L n +m 6L 

AtMo) ^ "W~ = 5L < 5L < 5L < 

rf T (r, di) _ d T (r, dp) + \d di\ 6L + |rf rfi| 

T(r ' l} ~ |rrfr| ~ |rdo| + |dodi| Sl + |cMi| 



On the other hand, 



A TM2 ) = ^>^^ = ^ = 1.4 
|rd 2 | \rd 2 \ 10L 



Hence, the delay of T is determined by d 2 : 

A(T) = max A T (r, t>) = A T (r, d 2 ) □ 

ti£S\{r} 

We define a special regular tree, the base tree To, which contains all regular edges except for the 
edges aid. Delay and cost of the base tree are as follows. 
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Lemma 4. The total edge length of the base tree is £{Tq) < 14. bL, and its delay is A(Tq) = 1.4. 



Proof. The total edge length and the delay of the base tree To are 

n 

£(T ) = \rd 1 \+Y,(\b l c t \) + \d 1 d 2 \ 

i=l 

1 ™ 1 " 

= 8L + - 2"* + 6£ < 14£ + - ^(7* + m) = 14.51 

i=l i=l 

a \ a ( r ] \ M1I + IM2I 8L + 6X 

A(T ) = A To M 2 ) = — ^ = = 1.4 □ 

Connecting d 2 to any point other than d\ will always produce trees with higher cost than the base 
tree: 

Lemma 5. If T is an irregular tree that contains an edge vd2 for v 7^ d\, then £(T) > £(Tq). 

Proof. Assume T contains an edge vd 2 , for v ^ d±, and consider the path Q in T connecting r and d±. 
There are two possible cases. First, assume that Q passes through d2- Then we have 

£{T) > £{Q) > \rd 2 \ + \d 2 dx\ = 10L + 6L = 16L > £(T ). 

In the second case, Q does not pass through d 2 . Since \vd 2 \ ^ jcio^l and £{Q) ^ \rd-\_\, we have 

£{T) > £{Q) + \vd 2 \ ^ \rdi\ + \d a d 2 \ = 8L + 3\/5L ~ 14.7082L > l(T Q ). □ 

We can now show that for 8 ^ 1.4, regular trees are better than irregular trees. 

Lemma 6. For every irregular tree T x with A(Ti) ^ 1.4, there exists a regular tree T* such that 
t{T x ) > i{T*) and A(T X ) > A(T*). 

Proof. If Ti includes an edge vd 2 with v ^ d\, then Lemma [5] implies the lemma with T* — Tq. We can 
therefore assume that T\ contains the edge d\d 2 , and no other edge incident to d 2 - 

For a spanning tree T of S, let 7r(T) denote the path connecting r and d\ in T. We define J(T) to 
be the set of indices i G {1, . . . , n} such that ir(T) contains all three vertices a,, &i, and Cj, but does not 
contain both edge ajCj and &jCj. 

Let T denote the set of spanning trees T of S such that d 2 is adjacent only to d\ in T, £(T) ^ £(Ti), 
and £(tt(T)) < ^(tt(Ti)). We have T x G T, so T ^ 0. 

We now pick a spanning tree T2 G T such that J(T 2 ) is minimal under inclusion. Let us suppose first 
that JiTi) 7^ 0. Then there is an i G J(T 2 ) such that 7r(T2) contains {ai,&i,Ci}, but does not contain 
both aiCi and 6^. Let it be the first of the three vertices encountered by tt(T 2 ), let v be the second one, 
and let w be the last one. One of the two edges incident to v in 7r(T 2 ) is different from a^Ci and biCi. 
Denote this edge by vp. Then \vp\ ^ 74 ^ \uw\. We obtain a new spanning tree T3 from T 2 by adding uw 
and removing vp. We have £(T 3 ) = £{T 2 ) + \uw\ - \vp\ £{T 2 ) s? £(Ti). Clearly %(T 3 )) < ^(tt(T 2 )), 
and so T3 G T. Furthermore, J(T^) C J(T 2 ) as i ^ ^(^a), a contradiction to the choice of T 2 . 

It follows that J{T 2 ) = 0. Let us define the set / c {1, . . . ,n} of indices i such that 7r(T 2 ) contains 
both edge ajCj and c,-6j. We define T* as the tree consisting of all regular edges, except that we remove 
aibi when i £ I, and that we remove a^c^ when i g" /. Fig. [2] shows an example of the regular tree we 
construct. Let Eq denote the set of edges a^Ci and Cibi for all i E I. By definition of /, Eq C 7r(T 2 ), and 
we define E\ — n(T 2 ) \ E . Let Eq and E{ be the projection of E and E\ on the y-axis. Then Eq U £7* 
must be equal to the segment rdi, and so £(Ei) ^ £(El) ^ — £(i?o)- ^ follows that 

£(tt(T*)) - \rdi \ — £{Eq) + £(E Q ) < £{E\) + £(E ) — £(tt(T 2 )) ^ £(tt(Ti)). 

This implies that A(T X ) ^ A Tl (Y,d 2 ) > A T .(r,d 2 ) = A(T*) by Lemma|3] 

For each point u G S \ {r}, let be the second vertex on the path from v to r in T 2 . We have 
^2) = E«es\M For the vertices on tt(T 2 ), we have J2v^(T 2 )\{r} \ v p(v)\ = £{ir(T 2 )). Since 
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Figure 2: Irregular tree T 2 (left) and its regular tree T* (right) 

^(^2) = 0, for each i $ I, one of the three vertices <Zj, bi, Ci is not on the path 7r(T 2 )- Since its nearest 
vertex has distance at least a*, we have 

^2) = E M y )l > W2)) + + l^il- 

»£S\{r) i0J 

On the other hand, we have 

i{T*)=t(-K{T*)) + Y,<Xi + \d 2 d 1 \ 

< £(tt(T 2 )) + J2 a i + Mil < i(T 2 ) < *(Ti). □ 

Correctness of the reduction. It remains to show that the constructed point set S has a spanning 
tree of small delay and small cost if and only if the original Knapsack instance had a positive answer. 

Lemma 7. The Knapsack instance has a positive answer if and only if there is a spanning tree T for S 
with delay A(T) < 1.4 + (W/1QL) and cost £{T) £(T ) - P. 

Proof. We first assume that the Knapsack instance has a positive answer. Let / C {1,2, ... , n} be a 
set of indices such that J2 i(£l Pi ^ P and Y^iei w t ^ W. Let T be the tree consisting of all regular edges, 
except that we exclude aibi for i £ I, and exclude a^Ci for i $ I. 
Then we have 

t(T) = £(T ) - ]T( 7l - ft) = £(T ) - < £ ( T o) - p ' 
iei iei 

d T {r, d 2 ) = \rd x \ + |dida| + + A - 7-0 = 14L + ^ < 14i + ^ 

iei is/ 

A(T) = A T (r,rf 2 ) < ~~Jq^~~ = 1-4+ (W/10L), 
and the claim follows. 

Assume now that T is a spanning tree for 5 with the given bounds. If A(T) < 1.4, then T must 
include an edge incident to d 2 other than d\d 2 and is not regular. But then Lemma [5] implies that 
£{T) > £(To), a contradiction. So A(T) ^ 1.4, and by Lemma [6] we can assume that T is regular. 

Since T is a spanning tree, it must include all regular edges, except that for each 1 ^ i ^ n, one of 
the three edges a^i, atCi, or biCi must be missing. We define / C {1, 2, . . . , n} to be the set of indices i 
such that T does not include the edge aibi. 







We have 

d T (r,d 2 ) = Ni| + |did 2 | + ^(a, + ft - 7i ) = 14L + ^ Wl . 

iei iei 

Since A(T) = A T {r,d 2 ) = d T (r,d 2 )/WL, we have 

Wi = d T (r, d 2 ) - 14L = 10XA(T) - 14L < 10X(1.4 + (W/WL)) - 14L = W. 



The cost of T is 



*(T) ^ *(T ) - £>i - ft) = t(T ) - 

iei iei 



and so 



£> > t(T ) - l(T) > £(T ) - (£(T ) -P) = P. 
iei 

It follows that the Knapsack instance has a positive answer. □ 

Reduction with integer coordinates. To complete our proof of Theorem [2] we need to construct 
a set of points with integer coordinates, such that the total number of bits is polynomial in the size 
of Knapsack instance. The construction given so far does not achieve this yet, since the points Cj 
are defined as the solution of a quadratic equation. We will therefore compute approximations c~i with 
|cj — Cj| < e, for an s to be determined later. The set of points obtained in that way will be denoted by 
S, which is the set of points r, dj, hi, c~i, do, d±, d 2 . In the following lemma, we bound by how much this 
approximation can change the tree cost and delay. 

Lemma 8. If T is a spanning tree on S and T is the corresponding tree on S, then \£(T) — £(T)\ < 12ns, 
and |A(T) - A(f)\ < 20ns. 

Proof. Let u, v be a pair of points in S, with it, v denoting the corresponding points in S. Since \uu\ < s 
and |w| < s, we have — < 2s. The tree T has 3n + 3 edges, and so \£{T) — £{T)\ < (6n + 6)e ^ 
12ns. 

Consider now X := dx{r,v), X := df(r,v), Y :— \rv\, and Y := \rv\. Since |w| < s, we have \Y— Y\ < 
s. The path from r to v in T passes through most n approximated points, and so \X — X\ < 2ns. Since 
the longest edge in T has length 10L and the path has at most 3(n + 1) edges, we have X < 30L(n + 1). 
We also have that Y ^ 4L by the construction of S. This means that X/Y < 7.5(n + 1). We get 

Y(X + 2ns) -X{Y -e) 
YY 

X(Y + s) -Y(X-2ne) 
YY 

And since Y ^ 1, 

2ns s X , 

+ ^ ■ — ^2ns + sx 7.5 (n + 1 < 20ns. □ 
Y y Y 

If the Knapsack instance has a positive answer, then S has a spanning tree T with A(T) ^ 1.4 + 
(W/WL) and £(T) ^ £(Tq) — P. On the other hand, if the instance has a negative answer, then this 
implies that for any subset of indices I C {1, 2, . . . , n} we have either J2iei Wi ^ W+lor ^ ieI Pi ^ P — 1. 
By Lemma[7j this means that any spanning tree T for S has either delay A(T) ^ 1.4 + ((W + 1)/10L) 
or cost £{T) ^ £{T ) -P + l. 

Let us set s = l/600nL. We approximate the c; with a precision of most e, resulting in the point 
set S. This set is the input to our problem, with a delay bound of 5 = 1.4 + (W/10L) + (1/20L), and a 
cost bound of K = £{f ) -P + 0.5. 

If the Knapsack instance has a positive answer, then by Lemma ^ there is a spanning tree T for S 
with A(T) s$ 1 A+(W/ 10 L) + 20ns = 1.4+(W/10£) + (1/30L) < 5 and l(f )-£(f ) > P-2Ans > P-0.5. 

On the other hand, if the Knapsack instance has a negative answer, then by Lemma[8j every spanning 
tree f for § has either A (f) ^ 1.4+ (W/10L) + (1/10L) - 20ns = 1.4 + (W/WL) + (1/10L) - (1/30L) = 
1.4 + (W/WL) + (1/15L) > S, or we have £(f ) - t(f) < P-1 + 24ns < P - 0.5. 
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In both cases, solving our single-source dilation-bounded minimum spanning tree problem correctly 
answers the Knapsack instance. 

By construction, the points r, a^, b i7 dj have integer coordinates. We construct the points Cj by solving 
a quadratic equation with an error of at most 2~ k , that is, with k bits after the binary point, where 
k is chosen such that 2 k > 600nL. Clearly k is polynomial in the input size. If we multiply all point 
coordinates in our construction and the cost bound by 2 k , then all points have integer coordinates. 
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